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In the past decade, our understanding of galaxy evolution has been revolutionized by the dis- 
covery that luminous, dusty, starburst galaxies were 1,000 times more abundant in the early 
Universe than at present. It has, however, been difficult to measure the complete redshift 
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distribution of these objects, especially at the highest redshifts (z > 4). Here we report a red- 
shift survey at a wavelength of three millimeters, targeting carbon monoxide line emission 
from the star-forming molecular gas in the direction of extraordinarily bright millimetre- 
wave-selected sources. High-resolution imaging demonstrates that these sources are strongly 
gravitationally lensed by foreground galaxies. We detect spectral lines in 23 out of 26 sources 
and multiple lines in 12 of those 23 sources, from which we obtain robust, unambiguous red- 
shifts. At least 10 of the sources are found to lie at z > 4, indicating that the fraction of dusty 
starburst galaxies at high redshifts is greater than previously thought. Models of lens geome- 
tries in the sample indicate that the background objects are ultra-luminous infrared galaxies, 
powered by extreme bursts of star formation. 

We constructed a catalog of high-redshift (z > 1) galaxy candidates from the first 1,300 
square degrees of the South Pole Telescope (SPT) 3 survey by selecting sources with dust-like 
spectral indexes in the 1.4 and 2.0 mm SPT bands 4 . A remarkable aspect of selecting sources 
based on their flux at millimetre wavelengths is the so-called "negative /^-correction' 55 , whereby 
cosmological dimming is compensated by the steeply rising dust spectrum as the source redshift 
increases. As a result, a millimetre- selected sample should draw from the redshift distribution of 
dusty starburst galaxies with little bias over the entire redshift range in which they are expected to 
exist. To isolate the high-redshift, dusty- spectrum galaxy population, sources with counterparts in 
the IRAS Faint Source Catalog^ (typically z < 0.03) were removed, and those with counterparts 
in the 843 MHz Sydney University Molonglo Sky SurvejP were removed to exclude sources with 
strong synchrotron emission (e.g. flat- spectrum radio quasars) that may have passed the spectral 
index cut. A sample of 47 sources with 1.4 mm flux density > 20 mJy and accurate positions 
were selected for high-resolution imaging with the Atacama Large Millimeter/sub-millimeter Ar- 
ray (ALMA). The ALMA spectroscopic observations targeted a sample of 26 sources, all but two 
of which are in the imaging sample (see Supplementary Information). These objects are among 
the brightest dusty- spectrum sources in the z > 0.1 extragalactic sky at millimetre wavelengths. 

Gravitationally lensed sources are expected to predominate in samples of the very brightest 
dusty galaxies because of the rarity of unlensed dusty, starburst galaxies at these flux level^^l. 
Massive elliptical galaxies, acting as lenses, will have Einstein radii as large as 2" and may magnify 
background galaxies by factors of 10 or more. To confirm the lensing hypothesis and determine 
magnifications, we imaged 47 SPT sources with ALMA at 870 fim in two array configurations, 
which provide angular resolutions of 1.5" and 0.5" (FWHM). A sample of these objects with IR 
imaging, spectroscopic redshifts, and resolved structure is shown in Figure 1. Integration times of 
only one minute per source are adequate to show that the sources are resolved into arcs or Einstein 
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rings - hallmarks of gravitational lensing - in most sources. For all sources for which we have 
IR and submillimetre imaging, as well as spectroscopic redshifts, the emission detected by ALMA 
coincides with massive foreground galaxies or galaxy groups/clusters, but is spatially distinct and 
at drastically different redshifts (see Figures 2 and S.3). Using a modeling procedure that treats 
the interferometer data in their native measurement space, rather than through reconstructed sky 
images, to simultaneously determine the source/lens configuration and correct for antenna-based 
phase error^, we are able to determine magnifications and derive intrinsic luminosities for our 
sources. Complete models of four lenses 11 , as well as preliminary models of eight more, indicate 
lensing magnifications between 4 and 22. After correcting for the magnification, these sources are 
extremely luminous - more than 1O 12 L and sometimes > 1O 13 L - implying star formation rates 
in excess of 500 M yr _1 . 

Obtaining spectroscopic redshifts for high-redshift, dusty starburst galaxies has been notori- 
ously difficult. To date, most spectroscopic redshift measurements have come from the rest-frame 
ultraviolet and optical wavebands after multi-wavelength counterpart identificatioiP^H. These 
observations are difficult owing to the extinction of the UV light by the dust itself, the cosmolog- 
ical dimming, and the ambiguity in the association of the dust emission with multiple sources of 
optical emission visible in deep observations. A much more direct method to determine redshifts 
of starburst galaxies, particularly at high redshift, is through observations of molecular emission 
associated with their dusty star forming regions. The millimetre and submillimetre transitions 
of molecular carbon monoxide (CO) and neutral carbon (C i) are well-suited for this purposed. 
These emission lines are a major source of cooling for the warm molecular gas fueling the star 
formation, and can thus be related unambiguously to the submillimetre continuum source 15 . Until 
recently, bandwidth and sensitivity limitations made this approach time-intensive. The combina- 
tion of ALMA - even with its restricted early science capabilities and only 16 antennas - and a 
unique sample of extraordinarily bright millimetre sources has changed this situation dramatically, 
allowing us to undertake a sensitive, systematic search for molecular and atomic lines across broad 
swaths of redshift space at z > 1. 

We conducted a redshift search in the 3 mm atmospheric transmission window with ALMA 
using five spectral tunings of the ALMA receivers to cover 84.2 — 114.9 GHz. For z > 1 at least 
one CO line will fall in this frequency range, except for a small redshift "desert" (1.74 < z < 
2.00). For sources at z > 3, multiple transitions (such as CO(J up > 3) and Ci( 3 Pi — ^ 3 P )) 
are redshifted into the observing band, allowing for an unambiguous redshift determination. We 
find one or more spectral features in 23 of 26 SPT-selected sources. The detections comprise 
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44 emission line features which we identify as redshifted emission from molecular transitions of 
12 CO, 13 CO, H 2 and H 2 + , and a Ci fine structure line. The spectra of all sources are shown in 
Figure 2. For 18 of the sources we are able to infer unique redshift solutions, either from ALMA 
data alone (12), or with the addition of data from the Very Large Telescope and/or the Atacama 
Pathfinder Experiment telescope (6). With the 10 z > 4 objects discovered here, we have more 
than doubled the number of spectroscopically confirmed, ultra-luminous galaxies discovered at 
z > 4 in millimetre/submillimetre surveys in the literature (of which just nine have been reported 
previousl^^E^H). Two sources are at z = 5.7, placing them among the most distant ultra- 
luminous starburst galaxies known. 

The SPT dusty galaxy redshift sample comprises 28 sources, as we include an additional two 
SPT sources with spectroscopic redshifts 22 that would have been included in the ALMA program 
had their redshifts not already been determined. Of the 26 ALMA targets, three lack a spectral 
line feature in the ALMA band. We tentatively and conservatively place these at z = 1.85, in the 
middle of the z = 1.74 — 2.00 redshift desert, though it is also possible that they are located at 
very high redshift or have anomalously faint CO lines. For the five sources for which only a single 
emission line is found, only two or three redshifts are possible (corresponding to two choices of 
CO transition) after excluding redshift choices for which the implied dust temperature, derived 
from our extensive millimetre/submillimetre photometric coverage (provided by ALMA 3 mm, 
SPT 2 & 1 mm, APEX/LABOCA 870 /im and Herschel/SPIRE 500, 350, 250 fim observation^ 1 ), 
is inconsistent with the range seen in other luminous galaxie^. For these sources we adopt the 
redshift corresponding to the dust temperature closest to the median dust temperature in the unam- 
biguous spectroscopic sample, as shown in Figure 2. 

The cumulative distribution function of all redshifts in this sample is shown in Figure 3. The 
median redshift of our full sample is z med = 3.5. The redshift distribution of SPT sources with 
mm spectroscopic redshifts is in sharp contrast to that of radio-identified starbursts with optical 
spectroscopic redshifts, which have a significantly lower median redshift of z med = 2.2, and for 
which only 15-20% of the population is expected to be at z > 3 2 . Part of this difference can 
be attributed to the high flux threshold of the original SPT selection, which effectively requires 
that the sources be gravitationally lensed. A much smaller total volume is lensed at z < 1 than at 
higher redshift, and, as expected, we do not find any such sources in the SPT samplePl However, 
if we only compare sources at z > 2 (the lowest confirmed spectroscopic redshift in the SPT 
sample), the median redshift of the radio-identified sample is still significantly lower (2.6) than 
the SPT sample, and the probability that both samples are drawn from the same distribution is 
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< 10~ 5 by the Kolmogorov-Smirnov test. A recently published survey 24 of millimetre-identified 
starbursts with optical counterparts determined from high resolution mm imaging and redshifts 
measured from optical spectroscopy or estimated from optical photometry found a median redshift 
of z med = 2.8. Again comparing the distribution of sources at z > 2, the probability that these 
objects and the SPT-selected sources are drawn from the same parent distribution is 0.43, indicating 
rough consistency between our secure redshift determinations and the distribution estimated from 
the optical methods. A full analysis of the molecular line detections, redshift determinations, 
residual selection effects, and a derivation of the intrinsic redshift distribution for the SPT sample 
is reported in a companion papei^. 

These 26 sources represent less than 25% of the recently completed SPT survey and catalog. 
This newly discovered population of high-redshift starbursts will enrich our understanding of ob- 
scured star formation in the early universe. Existing semi- analytic hierarchical models of galaxy 
evolutioip3251 have already had difficulties reconciling the number of z ~ 4 systems inferred from 
previous observational studie^^. The presence of two intensely starbursting galaxies at z = 5.7, 
1 Gyr after the Big Bang, in a sample of just 26 sources, demonstrates that significant reservoirs 
of dust and molecular gas had been assembled by the end of the epoch of cosmic reionization. As 
the millimetre-brightest high-redshift starbursts in the sky, the present sample will be key targets 
for ALMA studies of star formation physics at high redshift. The gravitational lensing of these 
sources provides access to diagnostic information from molecular lines that would otherwise take 
hundreds of times longer to observe, and effective source-plane resolution several times higher 
than can otherwise be achieved. 
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Figure 1: Figure 1: Near-infrared (NIR) and ALMA submillimeter- wavelength images of SPT targets. 
Images are S f/ xS f/ . We show 10 sources for which we have confirmed ALMA spectroscopic redshifts, deep 
NIR imaging, and well-resolved structure in the ALMA 870 iim imaging. The greyscale images are NIR 
exposures from the Hubble Space Telescope (co-added F160W and Fl 10W filters), the Very Large Telescope 
(K s ), or the Southern Astrophysical Research Telescope (K s ), and trace the starlight from the foreground 
lensing galaxy. The NIR images are shown with logarithmic stretch. The red contours are ALMA 870 fim 
imaging showing the background source structure, clearly indicative of strong lensing from galaxy-scale 
halos. In all cases, the contours start at 5 a and are equally spaced up to 90% of the peak significance, 
which ranges from 12 to 35. Spectroscopic redshifts of the background sources are shown in red in each 
panel. The ALMA exposures were approximately 2-minute integrations, roughly equally divided between 
the compact and extended array configurations. The resulting resolution is 0.5". SPT 0103—45 shows a rare 
lensing configuration of one lens and two background sources at different redshifts, one visible with ALMA 
and one with HST. SPT 0346—52, with a CO-derived redshift of z=5.656, is among the highest-redshift 
starbursts known. (See the Supplementary Information for more details.) 
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Figure 2: Figure 2: ALMA 3 mm spectra of 26 SPT sources. The vertical axis is observed flux density 
in units of mJy, with 30 mJy offsets between sources for clarity. Spectra are continuum-subtracted. The 
strong CO lines are indicative of dust-enshrouded active star formation. The spectra are labeled by source 
and redshift. Black labels indicate unambiguous redshifts (18), with the subset in bold (12) having been 
derived from the ALMA data alone. Five sources labeled in blue (5) are plotted at the most likely redshift 
of multiple options, based on the dust temperature derived from extensive far-infrared photometry. Three 
sources with no lines detected are placed at z = 1.85, in the middle of the redshift range for which we 
expect no strong lines, and labeled in red. Total integration times for each source were roughly ten minutes. 
The synthesized beam size ranges from 1" x 5" to 5" x 3" over the frequency range of the search, which is 
inadequate to spatially resolve the velocity structure of the lensed sources. 
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Figure 3: Figure 3: The cumulative redshift distribution of luminous, dusty starburst galaxies, as measured 
with different techniques. The SPT millimetre- selected sample, with redshifts directly determined from 
spectroscopic observations of the molecular gas in the galaxies, is shown in black. The existing samples 
of radio-identified starburst^^^^E^ with redshifts determined from rest- frame ultraviolet spectroscopy, 
are compiled in the blue distribution. The redshift distribution^ of millimetre-identified starburst galaxies 
in the COSMOS survey is shown in red/orange, though the majority of redshifts in this sample are derived 
from optical/IR photometry of the sources rather than spectroscopy, and therefore less certain. Sources at 
z < 1 were removed from the previous samples of starburst galaxies to better compare to the selection effect 
imposed on the SPT sample due to gravitational lensing. The distribution of redshifts for radio-identified 
sources is incompatible with the distribution for the sample presented in this work. This measurement 
demonstrates that the fraction of dusty starburst galaxies at high-redshift is greater than previously derived 
and that radio-identified samples were biased to lower redshift than the underlying population. 
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Sample Selection 

The dusty- spectrum sources targeted for the ALMA observations described here were found in the 
SPT survey. The full survey comprises 2540 deg 2 of mapped sky, but constraints of data analy- 
sis and followup limited the area available for target selection at the time of the ALMA Cycle 
deadline to 1300 deg 2 . The initial SPT target list was extracted from the SPT maps according to 
the procedure described in a previous paper 4 and the selection outlined in the manuscript. Before 
observing these sources with ALMA, we required that they have followup observations with the 
LABOCA 870 /xm camera on the Atacama Pathfinder Experiment telescope to improve the accu- 
racy of their positions. At the time of the ALMA Cycle proposal deadline, we had completed 
this followup for 76 sources within 1300 deg 2 of the survey area. We selected 47 sources for imag- 
ing and 26 sources for spectroscopy with ALMA; 24 of the 26 spectroscopic sources were also in 
the imaging sample. The sample selection targeted sources with the highest SPT 1.4 mm fluxes, 
subject to the restrictions of the ALMA call for proposals. The most important restriction was the 
requirement that sources be located within 15 degrees of each other on the sky. This should not 
affect the statistical properties of the sample, however, it merely prevented the observation of a 
complete set of SPT sources above a defined flux threshold. 

Detected line features 



The detected CO and C i line features are shown in Figure IS.4I Additional lines are detected in 
some spectra, including 13 CO transitions in two sources. However, the detection of both 12 CO and 
13 CO transitions in the same source does not break redshift degeneracies because both transitions 
are harmonically spaced; at a given frequency of detection, every pair of CO isotopic transitions 
of the same rotational level (J) will have the same observed spacing. Emission lines of H 2 and 
H 2 + are detected in the spectrum of SPT0346-52. 

ALMA imaging 

Continuum images from the ALMA observations at both wavelengths are shown in Figure |S3] for 
all 49 sources observed with ALMA in one or both of the 3 mm redshift search or the 870 /xm 
imaging projects. The positional coincidence between the bands confirms that the redshifts are 
derived for the same objects that are seen to show structures indicative of gravitational lensing. 
Nearly all sources are resolved at the 0.5" resolution of the 870 /urn data, most likely due to grav- 
itational lensing. Exceptions may be due to lensing by groups/clusters, with image counterparts 
that are either faint or too widely separated to be detected in the small (18") primary beam of the 
ALMA antennas at this wavelength, or because some of our objects have small image separations. 
Lensed dusty sources with similar image separations are already prominent in the literature, in- 



eluding APM 08279+525^3, which has three images separated by <0.4 // . Conclusive evidence of 
lensing in many objects, including the most compact, awaits IR imaging and spectroscopic red- 
shifts for the sources (and any candidate lens galaxies). The 23 objects for which we have the most 
complete data (the 26 sources of the 3 mm spectroscopic sample, less three without detected lines, 
less two without IR imaging, augmented by the two sources^ for which we had prior redshifts) are 
shown in Figure [ST6l 
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Figure S.4: The CO(3-2), CO(4-3), CO(5-6), CO(6-5) and Ci(l-O) emission lines observed with 
ALMA for 23 out of the 26 SPT sources, which were used to determine the source redshifts. 
The vertical axis is observed flux density, sources are offset from zero for clarity with the offsets 
specified in square brackets next to the source names. Redshifts marked with 'A' are confirmed 
using additional observations from other facilities, while redshifts marked with '?' are uncertain 
and are shown at the most likely redshift. SPT0452-50 has a single line, but is determined to be at 
z = 2.010 rather than z = 1.007 because the implied dust temperature for this source would be far 
lower than in any other source (13 K) were it at the lower redshiflpl. 
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Figure S.5: Continuum images of the 49 sources observed at 3 mm and/or 870 /xm wavelength 
with ALMA. The 47 sources for which 870 fim data were acquired are shown in greyscale, with 
the 3 mm images overlaid in red contours. Two sources from the redshift sample that lack 870 /im 
data appear as red contours on a blank background. Images are 10"x 10", the 870 /xm and 3 mm 
images have 0.5" and 5" resolution, respectively. The correspondence between the positions at the 
two wavelengths unambiguously links the lensing structure visible at 870 /xm to the 3 mm spectra. 
The 3 mm contours are plotted in units of 3a, starting at 3a for sources at S/N<15, and 5a for 
sources at S/N>15, except SPT 0457-49, where the contours are 3 and 4a. The grey scale stretch 
of each image is indicated in the lower left hand side of each panel and is roughly from —la to the 
peak value. 
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Figure S.6: Images of the full set of 23 sources for which we have ALMA 870 /xm 3 mm, or 
S ABOCA 350 jim imaging, deep NIR imaging, and a redshift for the background galaxy (including 
ambiguous redshifts). Except for SPT 0512-59, ALMA 870 /xm emission is represented with 5 red 
contours, spaced linearly from five times the image noise to 90% of the peak signal to noise, 
specified in the upper right of each panel. For SPT 0512-59, which lacks ALMA 870 /xm data, we 
show the ALMA 3 mm continuum contours. For SPT 2332-53, which lacks ALMA 870 /xm data, 
we show the APEX/SABOCA 350 /xm continuum contours. The redshift of the background source 
(zs) is specified in red. Greyscale images are near-infrared exposures from the Hubble Space 
Telescope (co-added F160W and F110W filters), the Very Large Telescope (K s ), the Southern 
Astrophysical Research Telescope (K s ), or the Spitzer Space Telescope (3.6 /xm) and trace the 
starlight from the foreground lensing galaxy. The images are shown with logarithmic stretch. 
When known, the redshift of the foreground galaxy (z L ) is specified in black. In nearly every 
case, there is a coincidence of the millimetre/submillimetre emission, determined by the redshift 
search data to arise at high redshift, with a lower redshift galaxy, a galaxy group, or a cluster. This 
is precisely the expectation for gravitationally lensed galaxies. Three cluster lenses are apparent, 
SPT 0550-53, SPT 0551-50, and SPT 2332-53, with two other systems lensed by compact groups 
(SPT 0113-46, SPT 2103-60). 



